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Abstract 



■ We study the radiation reaction on cosmic strings due to the emission of 

dilatonic, gravitational and axionic waves. After verifying the (on average) 

m 

J> ■ conservative nature of the time-symmetric self-interactions, we concentrate 

, on the finite radiation damping force associated with the half-retarded minus 

half-advanced "reactive" fields. We revisit a recent proposal of using a "local 

o . 

OO , back reaction approximation" for the reactive fields. Using dimensional con- 



tinuation as convenient technical tool, we find, contrary to previous claims, 
that this proposal leads to antidamping in the case of the axionic field, and 
^P' to zero (integrated) damping in the case of the gravitational field. One gets 

- ^ . normal positive damping only in the case of the dilatonic field. We propose to 

X 

$H ' use a suitably modified version of the local dilatonic radiation reaction as a 



substitute for the exact (non-local) gravitational radiation reaction. The in- 
corporation of such a local approximation to gravitational radiation reaction 
should allow one to complete, in a computationally non-intensive way, string 
network simulations and to give better estimates of the amount and spec- 
trum of gravitational radiation emitted by a cosmologically evolving network 
of massive strings. 
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I. INTRODUCTION 



Cosmic strings are predicted, within a wide class of elementary particle models, to form 
at phase transitions in the early universe IQ, 0. The creation of a network of cosmic 
strings can have important astrophysical consequence, notably for the formation of structure 
in the universe 0, @]. A network of cosmic strings might also be a copious source of 
the various fields or quanta to which they are coupled. Oscillating loops of cosmic string 
can generate observationally significant stochastic backgrounds of: gravitational waves 
massless Goldstone bosons p, light axions 0, [||, or light dilatons 0. The amount of 
radiation emitted by cosmic strings depend: (i) on the nature of the considered field; (ii) on 
the coupling parameter of this field to the string; (iii) on the dynamics of individual strings; 
and (iv) on the distribution function and cosmological evolution of the string network. It 
is important to note that the latter network distribution function in turn depends on the 
radiation properties of strings. Indeed, numerical simulations suggest that the characteristic 
size of the loops chopped off long strings at the epoch t will be on order of the smallest 
structures on the long strings, which is itself arguably determined by radiative back reaction 
10| , |[TT| . For instance, if one considers grand unified theory (GUT) scale strings, with 



tension fi ~ ^gvt^ gravitational radiation (possibly together with dilaton radiation which 
has a comparable magnitude ^) will be the dominant radiative mechanism, and will be 
characterized by the coupling parameter Gfi ~ (AcuT/^^pianck)^ ~ 10~^. It is then natural to 
expect that the same dimensionless parameter Gfi will control the radiative decay of the small 
scale structure (crinkles and kinks) on the horizon-sized strings, thereby determining also the 
characteristic size relative to the horizon of the small loops produced by the intersections of 
long strings: Aoops = « ct, with a ~ Fkink G/i, Fkink being some dimensionless measure of the 
network-averaged radiation efficiency of kinky strings [|T0|, [|r^, [|12|, [|1^]. If one considers 
"global" strings, i.e. strings formed when a global symmetry is broken at a mass scale fa, 
emission of the Goldstone boson associated to this symmetry breaking will be the dominant 
radiation damping mechanism and will be characterized by the dimensionless parameter 
fa / fJ'eSective ~ (}og{L / 6))~^ ~ 10^^, whcrc the effective tension /ieffective is renormalized by a 
large logarithm (see, e.g., 0). 

Present numerical simulations of string networks do not take into account the effect of 
radiative damping on the actual string motion. The above mentioned argument concluding 
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in the case of GUT strings to the hnk a ~ T^ink G/i between the loop size and radiative 
effects has been justified by Quashnock and Spergel ||TT| who studied the gravitational back 
reaction of a sample of cosmic string loops. However, their "exact", non-local approach to 
gravitational back reaction is numerically so demanding that there is little prospect to im- 
plementing it in full string network simulations. This lack of consideration of the dynamical 
effects of radiative damping is a major deficiency of string network simulations which leaves 
unanswered crucial questions such as: Is the string distribution function attracted to a solu- 
tion which "scales" with the horizon size down to the smallest structures? and What is the 
precise amount and spectrum of the gravitational (or axionic, in the case of global strings) 
radiation emitted by the combined distribution of small loops and long strings? 



Recently, Battye and Shellard |T^, [|T^] proposed a new, computationally much less inten- 
sive, approach to the radiative back reaction of (global) strings. They proposed a "local back 
reaction approximation" based on an analogy with the well-known Abraham-Lorentz-Dirac 
result for a self-interacting electron. Their approach assumes that the dominant contribution 
to the back reaction force density at a certain string point comes from string segments in the 
immediate vicinity of that point. They have endeavored to justify their approach by com- 
bining analytical results (concerning approximate expressions of the local, axionic radiative 
damping force) and numerical simulations (comparison between the effect of their local back 
reaction and a direct field-theory evolution of some global string solutions). 

In this paper, we revisit the problem of the back reaction of cosmic strings associated 
to the emission of gravitational, dilatonic and axionic fields, with particular emphasis on 
the "local back reaction approximation" of Battye and Shellard. Throughout this paper, we 
limit our scope to the self-interaction of Nambu strings, in absence of any non-trivial external 
fields. This problem can be (formally) treated by a standard perturbative approach, i.e. by 
expanding all quantities in powers of the gravitational^ coupling constant G. We work only to 
first-order in G. To this order, we first verify the fact (well-known to hold for self-interacting, 
electrically charged, point particles |jl6[) that the time-symmetric part of the self-interaction 



(i.e. the part mediated by the half-retarded 'plus half-advanced Green function) is, on the 



^ Because of our "gravitational normalization" of the kinetic terms, see Eq.(^.3|), the couplings of 
all the three considered fields are proportional to G. 
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average, conservative, i.e. that it does not (after integration) drain energy-momentum out 
of the string. As we are interested in radiation damping, this allows us to concentrate on 
the time-odd part of the self-interaction, mediated by the half-retarded minus half-advanced 
Green function. This "reactive" part of the self-interaction is (as in the case of point charges) 
finite. [By contrast, the time-symmetric self-interaction is (formally) ultraviolet divergent. 
This divergence is not of concern for us here because, as shown in Refs. [|l^], [|l^] and further 
discussed below, its infinite part is renormalizable, and, as said above, its finite part does 
not globally contribute to damping.] Contrary to the case of point charges, the reactive part 



is non-local, being given by an integral over the string. Following Battye and Shellard [|14 



15| we study the "local approximation" to this reaction effect. We find very convenient for 



this study to use the technique of dimensional continuation (well-known in quantum field 
theory). 

In the case of the axionic self-field, we find that the axionic reaction force defined by the 



"local back reaction approximation" of Ref. |]T4[ leads to antidamping rather than damping, 
as claimed in Refs. DTSf. We also investigate below the corresponding local approxi- 
mations to gravitational and dilatonic self-forces and find zero damping in the gravitational 
case, and a normal, positive damping for the dilatonic case. The physical origin of these 
paradoxical results is explained below (Sec. [IV F|) by tracing them to the modification of 
the field propagator implicitly entailed by the use of the local back reaction approximation. 
We show that, in the case of gauge fields, this modification messes up the very delicate 
sign compensations which ensure the positivity of the energy carried away by gauge fields. 
Thereby, one of the main results of the present work is to prove the untenability of applying 
a straightforward local back reaction approximation to gauge fields such as gravitational and 
axionic fields. However, this untenability does not necessarily apply to the case of non-gauge 
fields. Indeed, our work proves that the application of this local approximation to the dila- 
tonic field (which is not a gauge field) leads to the correct sign for damping effects. In this 
non-gauge field case, the argument (of Sec. [IV F| ) which showed the dangers of approximat- 
ing the field propagator for gauge fields, looses its strength. This leaves therefore open the 
question of whether the "local approximation" to dilatonic back reaction might define (de- 
spite its shortcomings discussed below) a phenomenologically acceptable approximation to 
the exact, non-local self- force. In this direction, we give several arguments, and strengthen 
them by some explicit numerical calculations, toward showing that the meaningful (positive- 
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damping) dilatonic local back reaction force can be used, after some modification, as a 
convenient effective substitute for the exact (non-local) gravitational back reaction force. 

This phenomenological proposal is somewhat of an expedient because it rests on an 
"approximation" whose validity domain is severely limited. However, pending the discovery 
of a better local proposal, we think that the incorporation of our proposed local reaction 



force (|5.1| ) should allow one to complete, in a computationally non-intensive way, string 
network simulations and to give better estimates of the amount and spectrum of gravitational 
radiation emitted by a cosmologically evolving network of massive strings. 

In the next section, we present our formalism for treating self-interactions of strings. We 



describe in Section |Tl| our results for the renormalizable, divergent self-action terms, and, in 
Section our results for the finite contributions to the "local" reaction force. In Section 
|V| we indicate how the local dilatonic damping force could be used in full-scale network 



simulations to simulate the dynamical effects of gravitational radiation. Section ^ contains 
our conclusions. Some technical details are relegated to the Appendix. 

As signs will play a crucial role below, let us emphasize that we use the "mostly positive" 
signature (— , +, +, +) for the space-time metric g^i, (/i, u = 0,1, 2, 3), and the corresponding 
(— , -|-) signature for the worldsheet metric 'jab (a, b = 0,1 being worldsheet indices). 

II. COSMIC STRINGS INTERACTING WITH GRAVITATIONAL, DILATONIC 

AND AXIONIC FIELDS 

We consider a closed Nambu string z'^{a"') (with = t, = a, < a < L) in- 
teracting with its own gravitational g^y{x^) = 7]^^, + h^v{x^), dilatonic ^{x), and axionic 
(Kalb-Ramond) B^^{x) fields. The action for the string coupled to gfj_^, ^ and B^y reads 

Ss = ~ [ Ai(v?) dA-^ [ dz^" A dz"" . (2.1) 
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Here dA = y/^d"^ a (with 7 = — det 7ab; 7^^ = g^u{z) da db z^ denoting the metric induced 
on the worldsheet) is the string area element and the dilaton dependence of the string tension 
\i can be taken to be exponential 

/i(¥^)=/ie2"^. (2.2) 



At the linearized approximation where we shall work the form (|2.2|) is equivalent to a linear 
coupling [i[}p) — fi{l + 2a!f). The dimensionless parameter a measures the strength of the 



coupling of if to cosmic strings (our notation agrees with the tensor-scalar notation of Ref. 
p9||), while the coupling strength of the axion field is measured by the parameter A with 
dimension (mass)^. Due to our "gravitational normalization" of the kinetic term of -B^jy, the 
link between A and the mass scale fa used in Refs. is 2G\^ = vr/^. 

The action for the fields is 



Sf = / d'x^ [7^ - 2 W - ^ e-^"^ H,,, H^^" 



(2.3) 



where H^^p = df^iB^^ + dyB^i^, + d^B^^, g = — det (g^u), and where we use the curvature 
conventions IZ^ypa = dpV^^ — ■ ■ ■ , T^nv = Tl^ ^pu- With this notation, a tree-level coupled 
fundamental string (of string theory) has a = 1 (in 4 dimensions) and A = /i. 

Everywhere in this paper, we shall assume the absence of external fields. More precisely, 
the background values of the fields we consider are g^^^^ = r]^^, (p^ = 0, B^^^ = 0. Our results 
are derived only for this case, by using (formal) perturbation theory around these trivial 
backgrounds. It is however understood, as usual, that one can later (e.g. for cosmological 
applications) reintroduce a coupling to external fields, varying on a scale much larger than 
the size of the string, by suitably covariantizing the final, trivial-background results derived 
here. Such an approximate treatment should be sufficient for the cosmological applications 
we have in mind. On the other hand, the methods used here are not appropriate for treating 
the general case of a string interacting with external gravitational and dilatonic fields of 
arbitrary strength and spacetime variability. To treat such a case, one would need a more 
general formalism, such as that of Ref. [^. Note, however, that the straightforward, non- 
explicitly covariant, perturbation approach to radiation damping effects used here is the 
string analog ^ of all the standard work done on the gravitational radiation damping of 



binary systems (see, e.g., |2T| for a review). 

In other words, our aim in this paper is to derive, consistently at the first order in the 
basic coupling constant G (see Eq. ( p.3|) ), both the fields generated by the string. 



^This analog is technically simpler because radiation damping appears at linear order for strings 
(which have a non-trivial, accelerated motion at zeroth order), while it is a non-linear phenomenon 
in gravitationally bound systems. 



6 



B^uix) = + GBl,ix) + OiG') , (2.4) 

and the (non-covariant) explicit form of the string equations of motion, written in a specified Q 
(class of) worldsheet gauge(s), 

/xr^,.(F-/-)=GJ-; + 0(G2). (2.5) 



The string action (2J.) can be written (using the Polyakov form) as 

5's = -|/ d^ae^'''^^r'daz''d,z''g^,-^ J d^ae'^^ d^z^ d,z^ B^, , (2.6) 

where the worldsheet metric ^ab must be independently varied and where = — 1, e^*^ = 1. 
The equation of motion of is the constraint that it be conformal to the induced metric 
lab = g^iiy{z) daZ^dbz'^. In the following, we shall often use the conformal gauge ^7 7"^ = 
^^^ab _ ^ab (^-y^^j-^grg ^^00 _ ^11 _ -)-i)^ j.e. t^^q shall choose the (t, cr) parametrization of 

the worldsheet so that 

i'^ g,, + z'^ z"'g^, = 0, z^ z'" g,, = . (2.7) 

Here z = doz = dz/dT and z' = diz = dz/da. Note also the expression, in this gauge, of 
the worldsheet volume density 

^ = g,,z'^z"' = -g^,z^z''. (2.8) 

Let us note that the string contribution to the energy-momentum tensor, 

T'^'^ = — ^ , (2.9) 
^/9 og^,u 



reads 

T^" = -^ [ d^ae'''^U'"'6\x- z(a)), (2.10) 

where / d'^x5'^{x) = 1 and 

U^" = -^Y^ da z" db z" = z^" - z'^ z'" (in conformal gauge) (2.11) 



We shall use the conformal gauge associated to the metric gi_iu{x) = rj^,^ + Gh]^^ + 0{G 
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is the "vertex operator" for the interaction of the string with the gravitational field (7^,^. [The 
vertex operators used here are the x-space counterparts of the fc-space vertex operators used 
in quantum string theory, e.g. U^^{k) = J (PaU^^{z{a)) exp(i /ea ^'^(o")).] The corresponding 
vertex operator for the interaction with the dilaton ip is simply the trace U = g^i, U^'^ , while 
the one corresponding to the axion B^^, is 

V^"' = -e"^ da db z"" = z"' - z" z'^" . (2.12) 

The exact equation of motion of the string can be written (in any worldsheet gauge) as 

= ^ = fig.^e^'^^daiVir'dbzn + , (2.13) 
where the quantity $^ is defined by 

$, = <i.;^ + <i.; + $J, (2.14) 

with 

<l>)^ = a e^"'^ U d^if - 2/i a e^"'^ g^^ f/"^ dp^ , (2. 15) 

$ J = ^ H^^p = ^ d^B^p + A d^Bp, . (2.17) 

Let us emphasize that, while and $^ are well defined, spacetime and worldsheet covariant 
objects, $1^, by contrast, is not a covariantly defined object (but the full combination 5Ss/5z^ 
of Eq. ( |2.13| ) is a covariant object). It can be noted that the sum of the dilatonic and 
gravitational contributions ( |2.15|) , ( p.l6| ) simplify if they are expressed in terms of the string 
metric (7*^, = e^°"-^ g^i, to which the string is directly coupled. Indeed, 

K + K = -/^ 9% K,[9%] U'^' = I f/"^ d,g:, - /i f/-^ d^g^,^ . (2.18) 

Except when otherwise specified, we shall henceforth work in the conformal gauge associated 
to the actual metric in which the string evolves (and not the conformal gauge associated to, 
say, a flat background metric 77^,^). In this gauge the equations of motion of the string read 

8^ = 0, with S^ = fig'^,r]''''dabz'' + <^^ = -fig'^A^''-z"") + '^f^- (2-19) 



When using this gauge, one must remember that the constraints ( p.7|) [which read the same 
when written in terms of the string metric g^^J\ involve the metric. These constraints read 
(in terms of the string metric) 

T,^, = 0, with T^^,^gl,{z)daz>^d,z^ ~]^r^ai>v'"'9'',Az)dcZ^d,z'' . (2.20) 

The constraints T^^ = are preserved by the (gauge-fixed) evolution ( p.l9|) . Indeed, it is 
easy to check the identity 

r^^' daT,% ^ d^z^ {g% r^^' da,z^ + + $ J) = d^z^ . (2.21) 

In the last step of Eq. ( |2.21| ) we used the algebraic identity dcZ^ = 0. When the gauge- 
fixed equations of motion are satisfied, i.e. when = 0, the constraints satisfy the conserva- 
tion law rj"''' daT^c = 0. This conservation law together with the algebraic identity r^"* T^f, = 
(i.e. Tqq = T^i), ensures that if T^^ vanishes on some initial slice r = tq , it will vanish 
everywhere on the worldsheet. This shows that the evolution equations ( p.l9| ) propagate 
only the physical, transverse degrees of freedom of the string. 

Up to this point, we have made no weak-field approximation. In the following, we shall 
limit ourselves to working with formal perturbative expansions of the form ( p. 41) , (|2.5|) . When 



doing this, it is convenient to rewrite the string equations of motion (|2.19|) in the explicit 
form 

£^ = -fiV^.Az''-z"n+^^.: (2.22) 
where the quantity JF^ is defined as 

= + (2-23) 

with 

vj/^ = {g^, e^^^- r^,,) (F - z'"') . (2.24) 

In the linearized approximation, the complementary contribution \E'^ to the equations of 
motion read 

= -/i (V + 2 « (/^ Vt^u) {r - z""') + Oi^C) . (2.25) 
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The total contribution JF^ to the exphcit (non-covariant) string equations of motion is not 
a covariantly defined object, it is a non-covariant, pseudo-force density. For the definition 
of a genuine, covariant force density see Ref. |2y], notably Eq. (41) there. To simplify the 



language, we shall however call, in this paper, the non-covariant combination a "force 
density" (in the same way that when doing explicit calculations of the perturbative equations 
of motion of binary systems it is convenient to refer to the right-hand side of the equations 
of motion "gravitational force".) 

Let us now write explicitly the weak-field approximation of the field equations deriving 
from the total action S"/ + Ss- Let us recall that we assume the absence of external fields, so 
that we work with perturbative expansions of the form ( p.4|) . When fixing the gauge freedom 
of the gravitational and axionic fields in the usual way (s'^'^r^^ = 0; B^^ = 0), the field 
equations derived from Sj read, at linearized order 

□^(x) = -An J (fa 5^{x - z{a)) + 0{G^) , (2.26) 

Uh^^{x) = -Atc J (fa Sj^ 5\x - z{a)) + 0{G^) , (2.27) 

UB^,{x) = -47r j (fa Sj, 6\x - z{a)) + 0{G^) , (2.28) 

where the corresponding linearized source terms are defined as 

E^ = aG^t/, s;, = 4G/it7^,, Sj, = 4GAr^,. (2.29) 

Here, U^y = U^i, — ^ rj^i, U, and, in the present approximation, the vertex operators entering 
these source terms are simply U^i, = z^z^ — z'^z^, U = rj^^" U^^, V^^ = z^z'^ — Zy z'^, where 
we freely use the fiat metric ri^^y to move indices. In the following, we shall consistently work 
to first order in G only, and shall most of the time omit, to save writing, the indication of 
the (9(G'^) error terms. 

The field equations (|2.26|) - (|2.28|) are classically solved by introducing the four dimensional 



retarded Green function 

act(x -y) = ^ ^^(x° - 2/°) 5((x - yf) ■ (2.30) 
UG,,,[x-y) = -b^'^\x-y). (2.31) 

This "retarded" Green function incorporates the physical boundary condition of the non- 
existence of preexisting radiation converging from infinity toward the string source. The 
(unphysical) time-reverse of Gret is the "advanced" Green function 
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1 



Gadv(a:-2/)= ^^(-(x°-y°))5((a:-yn 



(2.32) 

= G,,t{y-x). (2.33) 
Let us consider, as a general model for Eqs. ( p.26|) -( p.28D , the generic field equation 

UA{x) = -A-n j <fa^{a)5^{x- z{a)) . (2.34) 
Its most general classical solution reads 



+47r j da' dr' S(a', r') G,,t{x - z{a' , r')) + A,^,{x) 



(2.35) 



where At.(,^{x) is an "external" field, i.e. a generic homogeneous solution of the field equations 
(generated by far away sources). As said above, we assume in this work that Aext(a^) = 0. 
Applying the formula 5 {F{t')) = "^5 {r' — tq) /\d F{tq) / d tq\, where the sum runs over 

TO 

all the solutions tq of F{t') = 0, one can effectuate the integral over r' in Eq. (|2.35|) with 
the result 



Aret(x) 



da' 



[a ,r 



\n-z\ 



(2.36) 



Here, we have defined (x, a', t') = x'^ — z^ {a', r'), and Tret {x, a') as being the retarded (i.e. 
such that a;° — z'^ (ri.et(a;, a')) > 0) solution in r' of rj^yVL^^{T') Qy{T') = 0. In the following, 
we use also the quantity Aj-et which, after using the formula 

d^F dSiF) 



d^6{F{x,T'))=d,F6'{F{x,T')) 
and integrating by parts, can be written as 



{dF/d 



T' 



dr' 



di,A,et{x) 



da' 



1 d fQ^^{a',T') 



\n ■ z\ dr' 



n-z 



(2.37) 



f2.38) 



The corresponding results for the advanced fields are 

.,/S(a',r')\ 



^adv(3;) 



da' 



\Q-z\ 



(2.39) 



<9^^adv(a;) = J da' 



1 d fn^J:{a',T') 



\n ■ z\ dr' 



n ■ z 



(2.40) 



where Tadv{x^,(y') is the advanced solution of rif^i^Q^^^r') Q'^{t') = 0. Note that the scalar 
product Q ■ z is negative for t' = Tret and positive for r' = Tadv 
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III. PERTURBATIVE ON SHELL FINITENESS (AND RENORMALIZABILITY) 

OF THE STRING SELF-INTERACTIONS 



As said above, we consider the problem of a cosmic string interacting with its own, 
hnearized, gravitational, dilatonic and axionic fields. The equations of motion of such a 
string read 

/i77^, - = T^{z) + 0(6-2) ^ ^.^^^ jr^^^^ = dA'^\x%=, , (3.1) 

where the explicit expression of the linearized "force density" T^™ is a linear functional of 
the (linearized) retarded fields V5'^^'(a;), /i^'j^(x), -B]^°*(x), 

where 

^Minj^^ret] ^ ^ _ ^ ^^^^^^^ ^ 

^Blinj^^ret] ^ A ^^^^^^^ ^ ^ ^^^^^^^ ^ 

^'^[ip, h] = -/i (V + 2 « r^^.) (i;" - ;2"") . (3.3) 



The right-hand side of Eqs. (|3.3| ) is obtained by inserting the retarded fields A^et{.x)^ 
Eq. (p.36|) , and their first derivatives, (9^Ai.et(a;), Eq. (p.38|) , and by evaluating the result 



at a point = on the string worldsheet. The sources of the fields are given in terms 
of the string dynamics by Eqs. ( p. 291 ). Note that JF^ is a non-local functional of the string 



worldsheet whose support is the intersection of the worldsheet with the past light cone with 
vertex at the point . 

As in the case of a self-interacting point particle, the force J^fj,{x = z) is infinite because 
of the divergent contribution generated when the source point z^{T\a') coincides with the 
field point x^ = z'^{t, a). It was emphasized long ago by Dirac [jl6[, in the case of an electron 
moving in its own electromagnetic field, that this problem can be cured by renormalizing the 
mass, thereby absorbing the divergent part of the self-force. More precisely, Dirac introduced 
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a cut-off radius 6 around the electron and found a corresponding (ultraviolet divergent) self- 
force J-'^{6) = —{e'^/2 6)z'^ + -^R where JF^ is a finite (renormalized) contribution. If the 
mass of the electron plus its 5-surrounding depends on 6 according to 

^2 



■m{6) = mji — 



26 



(3.4) 



where denotes a finite, "renormalized" mass, the ultraviolet divergent equations of motion 
m{6) = give the finite result m-^'z^ = J-'j^. Note that the 5-dependence of m{6) (for 

a fixed mp,) is compatible with the idea that m{6) represents the total mass-energy of the 
particle plus that of the electromagnetic field contained within the radius 6 : 771(62)— m{6i) = 
+ J^^ d^x (Svr)^^ (e/r^y. Dirac also found that the remaining finite force was given by (using 
a proper-time normalization of r : = r]^,^ dz^ /dr dz"^ /dr = — 1) the sum of the external 
force J-'ltf and of a finite "reactive" self-force J-'H^^, 



'pfj- ■irA' _|_ •TTM -pi 

R — ext "T 



ext 



F'"'. z 



•'reac — o V-'^ret ^ sAw) '^v 



- e \z^\z-z)z 



(3.5) 



The analogous problem for self- interacting cosmic strings has been studied by Lund and 
Regge and Dabholkar and Quashnock |^ for the coupling to the axion field (see also 
), by Copeland, Haws and Hindmarsh p3 for the couplings to gravitational, dilatonic 



and axionic fields (and by Carter p5[ for the couplings to electromagnetic fields). There 
is, however, a subtlety in the calculation of the renormalization of the string equations of 
motion which led Ref. |^ (and us, in the first version of this work) to misinterpret their 
results, and propose incorrect values of the renormalizations of the string tension due to 
gravitational and dilatonic self-interactions. Our realization of this subtlety, was triggered 
by the work of Carter and Battye jl^, who were the first to get the correct renormalization 
of /i under self-gravitational effects, in 4-dimensions, by using a covariant approach to string 
dynamics 1^61 , |2^. We then obtained the correct renormalizations of [i under all three 
fields, and in an arbitrary spacetime dimension 0, by an effective action approach. The 
subtly which makes it delicate (but not impossible) to derive the correct renormalization of 



/i when working (as Ref. [24| and the present paper) directly with the equations of motion. 



Only the leading divergence was treated when n > 4. 
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at first order in G, and without adding external fields, is the following. In such a context, 
the perturbative string equations of motion ( |2.5D imply that z'^ — z^" is of order G, so that 
any first-order renormalization of the tension, /i = /i° + G /i^ + corresponds only 

to second order contributions {G {z^ — z'^) ~ G^/x^jF^ = 0{G^)) which are formally 
negligible at order 0{G) and, therefore, cannot be unambiguously read off such a first-order 
calculation. In other words, a first-order treatment without external fields can only prove 
that the string equations of motion are renormalizable by checking their on (perturbative) 
shell finiteness (i.e. the fact that all formally divergent first-order contributions vanish 
when using the zeroth-order string equations of motion i;^ — z^" = + 0{G)), but cannot, 
by themselves, unambiguously determine the renormalization of the string tension. For 
instance, the finding of Quashnock and Spergel |TT]] that the self-gravity effects vanish upon 
using the zeroth-order equations of motion to evaluate the first-order terms in Eq. ( |2.5| ) 
prove that they are renormalizable, but does not allow one to conclude that the self-gravity 
contribution to the tension renormalization, 6gfi vanishes. [It happens that 6gfi vanishes in 
4-dimensions fl^, but this vanishing is an "accident" which does not hold in other spacetime 
dimensions |]I8[ .] To be able to determine the value of the renormalization of /i one must go 
beyond a zero-background, first-order "on-shell" treatment of the string equations of motion. 
Essentially, one must work with a form of the string equations of motion which allow for the 
unambiguous introduction of an "external force" acting on the string. This is the case of the 



covariant-force formalism of Ref. [1 /l, as well as of the effective-action formalism of Ref. 



(where an extra force would mean an additional contribution to the total action.) In the 
present work, we do not really need the explicit value of the tension renormalization. We 
only need to check the renormalizability of the perturbative string equations of motion, i.e. 
the fact that all infinities vanish "on (zeroth-order) shell". To end up with clearer results, 
we shall, however, present a treatment in which the correct renormalization appear (because 
this treatment is action-based), and we shall renormalize them away by using, as external 
input, the results of Refs. [jl8[. 

Our starting point will be the explicit form, Eqs. (|0| )-(|0|), of the conformal-gauge, 
variational equations of motion SSg/Sz^ = 0. [We shall check below that the conformal- 
gauge constraints ( |2.7| ) (written with the full divergent metric) do not contain any divergent 
contributions (at linear order in G).] To give a meaning to Eqs. ( p.3|) when x —* z{T,a) 
we formally introduce an ultraviolet cutoff 6c in the cr'-integration giving the retarded fields 
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and their derivatives, i.e. we replace the integral over a full period of a', J^°^^ da', on the 
right-hand side of Eqs. (|2.36|) and ( p.38|) by J"'^" da' + XT+X^ da'. Later in this paper, we 
shall use a different way to introduce an ultraviolet cutoff, namely dimensional regularization. 
Dimensional regularization has the advantage of always keeping Lorentz invariance manifest. 
We have checked that both methods give the same results (see Appendix). In this section, 
we use the less sophisticated (5c-cutoff approach which allows a more direct comparison with 
other results in the literature. 

We then need the expansions in powers of a' — a and r' — r of all the quantities entering 
Eqs. (PI): 



1. / 



%{r', a') ^-{a'- a) < - (r' - r) - -ia' - af z'' - -ir' - rf I 



[a - a) (r - r) , 



■ z)(r', a') ~ - (r' - r) - \{r' - rf (i ■ ^) + ]^{a' - a)" {z" ■ z) 



+ [a'-a) (r'-r) {z ■ z') , 
Zf,{T', a') ~ i;, + (r -T)Zf, + {a' -a)z^. 



(3.6) 



(3.7) 
(3.8) 



At the order needed to extract the divergent part of the integrals ( p.36| ), ( |2.38|) (we shall use 
a more efficient tool below to extract the more complicated finite reactive part) it is enough 
to use 



T,,t{z^,a') = r - \a - a'\ + 0{\a - a'\') , 



(3.9) 



for the retarded solution of 77^,^ fi^(r') fi''(r') = 0. Inserting these results in Eqs. (|2.36|) 
(PI) we get 



Aret{z) = — y- log ^— ^ [2S] + finite terms. 



(3.10) 



d^Aj.ct{z) 



(i2)2 



+2S' z'^ - 2S 



-j:z, + ez';^ + 4Ez'^ 



finite terms . 



z"-z'' 



i2 



z ■ z 



(3.11) 



The rather complicated-looking terms proportional to {z' ■ z")/P and {z ■ z) / z^ in Eq. (|3.11| ) 
are, actually, "connection" terms linked to the fact that the source S is a worldsheet density 
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(conformal weight 2) rather than a worldsheet scalar (conformal weight 0). Let us associate 
to each source S a corresponding worldsheet scalar S, also denoted S, defined by 

1 



S = ^ 



(3.12) 



Here = (— det 7afe)^^^ is the area-density dA/(Pa, which reads, in conformal gauge: 
y/^ = z''^ = —P. One needs also to introduce the invariant ultraviolet cutoff S = 7^^^^ 5c = 
(z'^)^/^5c associated to the "coordinate cutoff" 6c- (In Sec. below and in the Appendix, 
we shall use a dimensional regularization method where the cutoff parameter e = 4 — n, and 
the renormalization scale Ar, are automatically Lorentz invariant). Then Eqs. ( p.lO|) , (|3.11 
simplify to 



log ( - ) [2 5] + finite terms . 



(3.13) 



Aret(^) = log [-Sz^ + Sz';^-2Sz^ + 2 S' z'^] + finite terms. (3.14) 



The result (|3.14|) for the regularized field derivative agrees with the results of Ref. 



as 



well as with the geometric prescription given in |^5[|. As a check on the above results one 
can verify that the divergent parts satisfy 

d 



OT 

d 

— Aret(^) = Z^" d^A,^,{z) . 

da 



(3.15) 
(3.16) 



To check these links one must use the following consequence of the conformal gauge con- 
straints = = g^v{z) da dhz" — }^ rjat Tj^^ g^iv{z) dc z^^ dd z" (written here in terms of the 
Einstein metric): 



= r{"^ d, = g,,{z) da v'' (d^d z" + da z") 



(3.17) 



As we assume everywhere in the paper a fiat gravitational background 77^1,, Eq. p.l7| ) implies 



z^{z, - z'') = 0{Kp) , z'^{z, - <) = 0{Kp) 



(3.18) 



so that tangential projections of z^ — z^" can be consistently neglected in first-order contri- 



butions such as Eqs. ( p.l5|) -( p.l7|) even if one is working "off-shell". 
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Because of the logarithmic divergence entering Eq. (|3.14|) we need to introduce, besides 
the invariant ultraviolet cutoff scale 6 (which can be thought of as the width of the cosmic 
string), an arbitrary, finite, renormalization length scale Ar. Then, we can define precisely 
the "infinite parts" (IP) of Aj.f.t{z) and dfj_ Aj.c,t{z) , i.e. the parts which blow up when 6-^0, 
by replacing in Eqs. (|3.13|) , (|3.14|) the logarithm by log (Ar/(5), and by discarding any other 
finite contribution. To apply this definition to the three fields ip, h^y, and B^y, we need to 



use the corresponding sources, Eq. (|2.29| ). For instance, we have 



^nKli^)) = ^2S;, log (^) = 8G/i-Lt7,. log ) , (3.19) 



IP(^^^*(^)) = -4aG/i log . (3.20) 

Using the easily verified identities satisfied by the vertex operators, 

U^yU^" ^U^yU^" --U^ = [/^^[/- = 0, U^,U =0, (3.21) 
V^yV^" = -2{z^f , V^,V^'' = -2-i'f + 2fz'^-A{i-zf-A{z'-zf, (3.22) 

we first see easily that the divergent contributions, YP{h^jj) (i^ z'^ + z^' z'^') and IP{h^y) z'^', 
to the constraints (p.7|) vanish. The use of the identities ( p.21|) , (|3.22| ) allows also to simplify 
the expression of the terms linear in the field derivatives entering Eqs. ( |3.3| ). We obtain 



IF (8,^) = -2a Gfi{z,-z;)- log [^^j , (3.23) 

IP([/"^ d^hp,) = -SG^iiz, - zD log , (3.24) 

IP(f/"^9^M = 0, (3.25) 

W{V^P d.B^p) = 8GX (z, - z;) log , (3.26) 

IP(y"^ d^Bfs,) = -80 X (z^ - z;) log . (3.27) 

We have now in hands all the results needed to derive the infinite contributions to the 
right-hand side of the string equations of motion (|3.1|) . More precisely, one obtains for each 
separate contribution in Eq. (|3.2| ) 

IP($):) = -Aa' G /x^ {z, - <) log (^) , (3.28) 

IP($;) = 8G /i^ {z, - <) log (^) , (3.29) 
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IP(<|.^) = -AG X' {z, - <) log ( ^ ) , (3.30) 



t 

Ar 



\nK) = ^a'G^^\i,-zl)\og(^] , [3.31] 



IF{^';) = -8G/i^(i, - z^) log j . (3.32) 
Adding up all the terms leads to 

ip(j-,) = c(i;,-^;)iog(^) , (3.33) 

with 

C = +4a^Gfi^ -4GX\ (3.34) 



The crucial point in the result ( |3.33D is that the divergent contribution to the equations of 



motion is proportional to the zeroth-order equations of motion. In our present perturbative 
treatment the formally infinite contribution ( p.33 ) is of second order in G and can be ignored. 



As we said above, this property of perturbative on shell finiteness of the equations of motion 
proves their renormalizability but cannot, by itself, determine the physically correct value of 
the renormalization of fi. At this point, we can, however, use the results of Ref. [|1^], where 
we showed that the "bare" (regularized but not renormalized) string tension fi{6) appearing 
in the original ultraviolet-divergent action must depend on the UV cutoff 6 according to 

Ar 



fi{6) =fiR + C log j , (3.35) 

where is the finite, renormalized tension, and where the ("beta function") coefficient C 
is precisely given by Eq. ( p.34[ ). [C contains only contributions coming from dilatonic, G^ = 



+4a'^ G fi"^, and axionic, Gb = — 4G'A^, self-interactions. The gravitational contribution 
vanishes (in 4-dimensions) 



Let us define the "renormalized" value of any "bare" (i.e. cutoff-dependent), logarithmi- 
cally divergent, quantity Q{S) as its "finite part" (FP), i.e. the difference between Q{6) and 
its "infinite part" (defined above as the term oc log(A/j/5) in Q{S) ) 

= FP {Q{6)) ^ hm [Q{6) - IP {Q{6))] . (3.36) 

0— »U 

Using this definition, and formally inserting Eq. ( |3.35|) into the bare equations of motion 
( p.l|) , namely 



+ iz" - z^") =[^iR + C log 



ri,. {i" - zn 



C {-z, - z'l) log ( ^ ) + + , (3.37) 



J-^(5) = IP(J-^(5))+^; 

'A 
~b 

we see that the terms proportional to log coming from the renormalization of and 
those coming from the renormalization of T^{S) are identical (even if we were working off 
shell), so that the equations of motion can be rewritten in the renormalized form 

/in - zn = + 0{G') . (3.38) 

This simplification between the same C log(A/j/(5) {z^ — z'^ contributions on both sides of 
the equations of motion is due to the fact that we have been working with the direct, Euler- 
Lagrange variational equations 6Ss/Sz^, i.e. with a form of the equations of motion which 
is ready to receive an additional "external force" 6S'/6z'^, as the variational derivative of 
an additional piece S' in the action. Had we worked with another form of the equations of 
motion, say 

gf" = = g^"' £^ = -fi{5) (F - z^") + $^ , (3.39) 

oz'^ 

with = g^^" 6^"^°"-^ the infinite part of the linearized field-contribution $^'™[9^Arct] 
would have been identical to rj^^ ^J^'^fc^^Aret], with given in Eq. (|3.3|) above. In such a case, 
Eqs. ( p.2^ )- (|3.30| ), show that the infinite part of would not have matched the infinite 
contribution IP(/x(5)) {z^ — z^") = C \og{/S.R/5) {z^ — z'^"). This apparent discrepancy is, 
however, not at all a sign of inconsistency of the type of non-covariant perturbative equations 
of motion we have been using. Either one works on shell, and all the formally infinite 
terms can be consistently neglected as being of order G^, or one introduces an additional 
mechanical interaction of the string, e.g. through the addition of a new piece S'{z^, . . .) in 
the action, in which case the zeroth-order string "mass shell" is modified, and we must take 
into account the new infinite terms coming from the extra contribution g^'^ g-scn^ 6S'/6z'^, in 
which g^"" e'^"^ = 7]^"" - lP{h^''') - 2a r/^^ IP{<^) + finite. 

Let us finally note that the logarithmic renormalizations ( 3.35|) , ( |3.36 ) introduce a de- 



pendence of the renormalized quantities upon an arbitrary, renormalization length scale A/j. 
[By definition, the bare (regularized) quantities fi{S), Q{6), do not depend on the choice of 
A/j.] For instance, we see from Eqs. (|3.33| ), ( p.35| ) that 
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M^'^) = M^r) - C log (^^^ , (3.40) 

Ki^K) = -^.^(Ar) -C{z,- z;) log (^^) . (3.41) 

It is however, easily seen that the content of the renormalized equations of motion ( p.38| ) is 
left invariant (at first order in the field couplings) under a change of Ar . [This invariance still 
holds in presence of an additional (finite) contribution 6S'/6z^ to the equations of motion.] 
As we work only to first order in the field couplings, note that the quantity fi appearing 
in C, Eq. ( |3.35|) , can formally be considered as being a renormalized value, rather than the 
bare one, thereby leading to the renormalization group equation d fi^/d log Ar = — C(/iR). 
[The non-renormalizability of the gravitational interaction makes it delicate to extend this 
argument to higher orders in G. By contrast, if we consider only a canonically normalized 



axionic field, with couphng VGA = ^11/2 fa, C does not depend on /i and the first-order 
renormalization result is exact.] 

Finally, we note that the axionic contribution Cb = — 4G to C agrees with the result of 
previous dynamical calculations |^4|, [|ll|], [Q and ||15|, while the dilatonic contribution 



= +4a G fx disagrees with Ref. p4| which proposed a vanishing dilatonic contribution 



IV. RENORMALIZED FORCE DENSITY AND THE LOCAL BACK-REACTION 

APPROXIMATION 

A. Renormalized equations of motion 

In the previous section we have shown that the perturbative equations of motion (in 
absence of external fields) could be written, at first order in G, in the renormalized form 

mVf.. {'^ - = Tf^idA^ + O(G^) , (4.1) 
where the R.H.S. is the sum of three renormalized contributions 

j'f'^ldA,,] = j-;^ii-[av7.et] + .F;^i^°[9/i,et] + .F^^^""[ai?,et] , (4.2) 

with 
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_^hRiin(^^) = I U-P {d.Kp)'' - U^^ {d^hp.f , (4.4) 



^BRiin^^^) = ^ V''^ {d.B^py" + A (d^Bp.f . (4.5) 



Here (dfj^A^et)^, with Aret = (v^''*'*, ^JTi^, -B^it*), denotes, as defined by Eq. ( |2.36| ), the finite part 



of the logarithmically divergent retarded integral ( |2.38| ). Note that, due to the absence of 
external fields, the supplementary contribution \E'^ to JF^, in Eq. ( |3.3| ), is negligible, being of 
order because hf^^, + 2a(p = 0{G) and {z^ — z'^) = 0{G). [Both the infinite part and the 
finite part of are 0{G^).] 

The expressions ([4.3[)-(^75D are linear (non-local) functions of the field derivatives. Fol- 



lowing Dirac [IC] it is useful to decompose any field Aj-^tix) in two parts: 



Axet{x) — Agjjj^i^x) + Aj-eac 

(x) , (4.6) 

A,y^{x) = i (Aret(x) + Aadv(a;)) , (4.7) 
^eac(a;) = ^ A,^a{x) = ^ {A,,t{x) - Aadv(a;)) . (4.8) 

Note the definition of two fields, Areac and Aj-ad? differing by a factor 2, associated to the 
difference A^^t — ^adv Both fields play a special role in the discussion below. They are both 
finite, as well as their derivatives, when considered at a point x = 2; of the source. Therefore 
the contribution to the self-force corresponding to A^^ac is finite and does not need to be 
renormalized. Hence, we shall dispense in the following with the label R when considering 
J-'^{Ay.eac) ■ To simplify the notation we henceforth drop the label "lin" on JF^, and freely 
move indices by rj^,^ because we shall consistently work only to first order in G. [As said 
above, in the present (first-order, no-external-field) approximation, we could even formally 
dispense with renormalizing jF'^(Asym) because the divergent contributions (|3.28| )- (|3.3CI|) are 
0{G{z^ — z'P) = Oi^G"^). But, for clarity we continue to work with jFj^(Asym)-] 

B. Reactive part of the self- force 

Let us first prove why, very generally, in the decomposition of the force corresponding to 

Eq. 
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•^R(^ret) — •^R(^sym) + •^'^(^reac) = -^Rsym + -^reac > (^-Q) 

the term ^/^ac considered as defining the full radiation reaction force, responsible for 

draining out of the mechanical system on which it acts (the string in our case) the energy 
lost to infinity in the form of waves of the A field. Indeed, for any field (in the linear 
approximation) we can define a field (pseudo) energy-momentum tensor T^'^{A) which is 
quadratic in (the derivatives of) A. The total energy tensor T^'^ — T^^ + ^ (where Tj^'^ 
denotes the energy tensor of source including any possible field-interaction energy localized 
on the source) is conserved: Q — di, T^". This leads to the equations of the source: d^, T^^ — 
Fi^{A) where F'^{A) = —d^, T^'^{A) represents the spacetime (rather than worldsheet) version 
of the force density acting on the source. [We work here with the bare force density.] Let 
us consider, as a formal simplification, the case where the coupling between the source and 
the field A is (adiabatically) turned off in the far past and the far future. [This means, in 
particular, that any possible field-interaction energy localized on the source vanishes in the 
far past and the far future.] Then the energy-momentum lost by the source during the entire 
interaction with the field, -Pfiogt — — J d^x [T^^{+oo) — Tj^°(— oo)], can be written as 

n^iost = -Jd'x 9^ = -jd'x F^{A) = P;^,,,,, , (4.10) 

where -Pjgained — + / [^/'^(+°°) ~^/°(~°°)] is energy-momentum gained by the field. 
When applying this result to the usual interaction force F^(A,-et) — —d^Tf^ (^ret) one has 
zero energy in Aret in the far past, so that P^^iost = / cPxT^^ {A^et{t = +00)). The field 
energy momentum tensor Tj^{A) is quadratic in the field and can always be written as the 
diagonal value of a symmetric quadratic form T'^^^A) = Qf^^^AjA). It is easy to see that 
the generic structure F'^{A) = —di,Tf^{A) — Sa ■ d A, where Sa is a source term for the 
field A, and where the dot product denotes some contraction of indices, is generalized, when 
considering Q'"'' to: -d^ (5''''(^i, A2) = \ [Sa^ • 9 A2 + • d A^\. We can apply this to the 
case where Ai — Aj-^^ and A2 — Aj.^ — ^^et ~ ^adv (for which Sa2 — 0) with the result: 

]^S-dA^^S- a Aeac = F^{A,^) = -d, g'^^(Aet, Aad) , (4.11) 

where S is the usual source, and F^{Aj-ca,c) the result of replacing Ajct by Aj-eac = \ ^rad in 
the usual force density. Integrating the latter formula over spacetime gives 

- / d'xF^{A,,^) = / d3x[g'^°(A,et,^ad)| -Q'^°(^et,^ad)| ]■ (4-12) 

J J \t=-\-oo \t= — oo 
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Again, one has zero energy from the far past contribution (because ylret(— cxd) = 0), while 
the far future contribution is simply, thanks to A^^vi+oo) = 0, Q^^{A,.ct, Aj-d) = T'^^{Aj.ct) 
so that 

- J d^xF'^iA,^) = J d'xT^\A,,i+oo)) = -J d^xFf^iA,,) = Pi'i,,^ ■ (4.13) 

This proves, for any field treated in the linear approximation, that the contribution to the 
self-force due to Areac contains, when integrated over time, the full effect of radiation damp- 
ing, ensuring conservation with the energy-momentum lost to radiation. The contribution 
F'^(Ai.eac) can be called the "reactive" part of the self-force F^{Aj.et)- 

Summarizing the results at this point, the renormalized self-interaction force (returning 
now to the worldsheet distributed force density) can be written as 

jr^ = FFJ^^ = FFT!^^J6) + T^,^ , (4.14) 

where FP denotes Hadamard's Finite Part ("Partie Finie") operation (i.e., in our case, 
the result of subtracting a term oc log(5/AR) from the ultraviolet-cutoff integral JF^(5) = 
Sai^^" do'[. . ] + ^^]^^ do'[. . .]). Note that only the symmetric contribution, obtained by 
replacing Aj-^t by A^yra = \ (^ret + ^adv) in the force density, needs to be renormalized (and, 
as we said above, one can even formally dispense with considering this renormalization). 
This symmetric contribution does not contribute, after integration over time, to the overall 
damping of the source. The finite reactive contribution J^^^^ = ^^(^reac) embodies (on the 
average) the full effect of radiation damping. 

The advantage of the above decomposition is to isolate, very cleanly, the radiation damp- 
ing force from the other non-cumulative, self-interactions. Its disadvantage is to write 
the non-local, but causal self-force FP J-''^{A^ct) as a sum of two acausal (meaning future- 
dependent) contributions. Indeed, both FP(jFg'^j^) and J^^ac given by integrals whose 
support is the intersection of the worldsheet with the two-sided light cone with vertex lo- 
cated at z^. In principle one can work directly with the full, causal (as done, e.g., in 
Ref. [jlll), but this is computationally very intensive. [A simplification, used by the lat- 
ter authors, and mentioned above, is that the self- force JF^(5) becomes, as is clear from 
Eqs. ( |3.28| )-( p.30| ), finite as 5 ^ when evaluated on free-string trajectories, satisfying 
— z^" = 0.] We shall follow Refs. [Q, [|1^ in working only with the (finite) reactive force 
jF^'^ac in trying to define a simple local approximation for it. 
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C. Local back-reaction terms in dimensional regularization 



The reaction force J-'^^a.c is linear in d^j, A,.eac{z), which is itself given by the following 
integral 



(4.15) 



with 



1 d fna^{a',T' 



\ri ■ z\ dr' \ n ■ z 



1 d (n^E{(T',r') 



|f2 ■ il dr' \ ■ z 



'''adv 



(4.16) 



The integrand B^{a') is the finite difference between two terms that blow up when a' ^ a 
{a being such that z = z(ct, r)). When a' is well away from a (say, for long, horizon-sized 
strings) B^{a') is expected to decrease roughly as the inverse spatial distance |f2 ■ i|, i.e. 
roughly as \a' — a\^^. In other words, a very rough representation of the typical behaviour of 
B^{a) is B{a') ~ (2(a' - a))-^ [/(r - {a' - a)) - /(r + {a' - a)], where the "effective source 
function" /(r) is expected to oscillate as r varies. If we think in terms of one Fourier mode, 
say /(r) = /^^e"*'^'^, these considerations suggest that the field derivative d A is roughly 
given by an integral of the form 



dA= da'B{a') r^if^e 



hoo 

' da' 

-oo 



smo; a 



o- 



(4.17) 

[a' - a) 

The latter integral is equal to vr, so that one can finally replace the oscillatory and decreasing 
integrand B{a') by an effective (5-function, Bg^^a') = B{0) A 6{(t' — a), with (in our example) 
B{0) = — /(O) = iujfi^ and A = tt/uj, or, in other words, dA = J da' B{a') is replaced by 
A 5(0). The analogous proposal of replacing the complicated, non-local integral ( [4.15| ) giving 
Aj-eac simply by the local expression 



[d^^A, 



1 local 



A 5^(0) 



(4.18) 



where A is some length scale linked to the wavelength of the main Fourier component of the 
radiation, was made by Battye and Shellard [0], [|l5l (see also |2^). In effect, this proposal 
is equivalent to replacing the a'-extended source S(r',cr') by the a'-local effective source 
AS(r',cr) 6{a' — a). One of the main aims of the present paper is to study critically the 
consequences of this proposal. 
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Though this "local back reaction approximation" drastically simplifies the evaluation 
of the reaction force ^^ac there remains the non-trivial analytical task of computing the 
a' — s> limit of the difference between the two complicated (and divergent) terms making 
up B^{a'). We found very helpful in this respect to use dimensional regularization, i.e. 
to use, instead of the normal (singular) four dimensional Green's functions ( p.30|) , (|2.32 ), 
their analytic continuation to a spacetime of (formal) dimension n = 4 — e. [We shall 
keep computing the index algebra in 4-dimensions. This is allowed here because our use of 
dimensional regularization is, simply, a technical trick for computing the finite object B^{a).] 
This technique is well known to be quite useful in quantum field theory, but it (or, at least, a 
variant of it) has also been shown long ago to be technically very convenient in the classical 
theory of point particles [g§, 0, g, . 



Riesz I^S] has shown that the retarded and advanced Green's functions in dimension 
n = 4 — e read 

G^:l {x~y) = -1- {-{x - yff-^ e{~{x - yf) 9i±{x' - y°)) , (4.19) 

adv \ / 

with Hn{2) = 2 7r'^ r (^) and 

□ G(")(x -y) = -5"(x - y) . (4.20) 

Note that, when e = 4 — n 0, the coefficient appearing in Eq. ( [4. 191) becomes 

e 



To save writing, we shall neglect in the following the factor 1 + (9(e) in Eq. ([4.21| ) which 
plays no role in the terms we consider. Then, we write the retarded solution of our model 
field equation (2.1) in dimension n = 4 — e as 

A,,t{x) =e I da' dr' S (1^')^ e{n^) , (4.22) 

J J — CO 

where fl^ = — (a; — z(r', cr'))^. (Note the inclusion of a minus sign so that Q'^ > within the 
light cone). Again neglecting a factor 1 + 0{e), the field derivative reads 

/PTrct 
da' c/r'Sfi^(^]2)-"/2^(^]2). (4.23) 
J —oo 

Using some efficient tools of dimensional regularization (which are explained in Appendix 
A) we get our main technical results: the explicit expressions of the reactive field, and its 
derivatives, in the local back reaction approximation 
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[^reac(^)] 



local 



local 



A 

i2 



A 



S - S 
1 



z ■ z 



2;^ 



2^2 



z^ 



(4.24) 



(4.25) 



z^ J 6 \ ) \ Z^ 

Some simplifications occur if we introduce, instead of tlie worldsfieet density S, tlie corre- 
sponding worldslieet scalar 5 = 2/ ^7. We find 



-S-S 



z ■ z 



[9^A,eac(^)]'°''' = ^ Is 



Z ■ Z 



2z 



z ■ z 



2; ■ z 



iS z^ . 



'Z^ o 



(4.26) 



(4.27) 



Note that Eqs. (|42^ ), ( ^^25] ) and Eqs. ( |42^ ), (|]2^ satisfy the compatibility condition 
d^A = A, but because of the lack of worldsheet covariance (broken by the introduction 
of A) the analog condition for z' is not verified. 



D. Dilaton radiation reaction 



Let us first apply our results to the case of the dilaton field (f, which has not been 
previously studied in the literature. The corresponding (worldsheet scalar) source is then 
simply 



= —T.^ = —aGfiU = -2aGfj. . 
Sip being a constant, the preceding formulas simplify very much: 



local 



2aG/iA 



z ■ z 



i2 



[5/.V5reac(2:)] 



local 



A 



1 ... 2 z-z\ 

- Zn -\ Zn TT- — 2Z 



z ■ z 



/II .9 



+ Zn 



z ■ z 



Inserting these results in the dilaton self-force ( ^4.3| ), we get 



local = _c,2G'/i2 A 



z ■ z 



z ^ z^ 



.2 



z ■ z 



(4.28) 

(4.29) 
(4.30) 

(4.31) 



z^ J \ z^ ^ 

For notational simplicity, we henceforth drop the label "local" on the local approximations 
to the reactive forces. Consistently with our choice of conformal gauge (which, in the case 
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of the dilaton coupling, is the same as in flat space, see Eq. ( p.l7|) ), we see that the reaction 
force ( |4.31 ) is orthogonal to the two worldsheet tangent vectors, and z'^: 



Tl = ^ = z'^" T^. (4.32) 

Let us now show that the putative, local approximation to the dilaton reaction force, Eq. 
( [4.31 ) conveys some of the correct physical characteristics expected from a radiation damping 



force. In particular, let us check that the overall sign of Eq. ( |4.31| ) is the correct one. First 
we remark that we can work iteratively and therefore consider that the reaction force ([4.31|) 



and its integrated effects, can be evaluated on a free string trajectory. In other words, 
when evaluating the total four momentum lost, P'f^^ = P]™(— oo) — P]™(+oo) (with, say, 
the convenient definition -P!j^(r) = Jq^ rfcr /iR i;^(cr, r)), by the string under the action of 

h I tB 



PlT = -l dadrT,, (4.33) 

we can insert a free string trajectory on the right-hand side of ( [4.33| )p]. This being the case, 
we can now further restrict the worldsheet gauge by choosing a temporal conformal gauge, 
i.e. such that t = -2;°(r, cr) = r. Geometrically, this means that the r = const, sections of the 
worldsheet coincide with x° = const, space-time coordinate planes. [The choice = r is 
consistent for free string trajectories because for them z'^ — z^" = 0.] In this gauge, we have 

i° = 1 , -i^ = 1 _ = z'^ i . z = V ■ V z' -z = z' (4.34) 

where we have introduced the 3-velocity v = z. The zero component of Eq. ( [4.33| ) then 
reads 



J^^ = +^a'Gf^'A^^. (4.35) 



Assuming that the scale A is constant, we can integrate by parts and write for the total 
energy lost by the string 



^Strictly speaking the integral in Eq. (4.33) is infinite because free string trajectories are periodic. 
The meaning of Eq. ( 4.33| ), and similar integrals below is to give, after division by the total 
coordinate time span r, the time-averaged energy-momentum loss. 



27 



3 



4 



V2 



+ 2- 



V ■ V 



(4.36) 



_1-V2 (1-V2)2_ 

The integrand of Eq. ( [4.36 ) is positive definite, ensuring that the reaction force ( 4.31|) has 
the correct sign for representing a radiation damping force. 

We can further check that the total 4-momentum lost by the string is, as it should, 
time-like. First, let us note that the relation 



(4.37) 



shows that, as far as its integrated effects are concerned, the dilaton reaction force ( p.l5|) is 
equivalent to 



V-'reac ^ ^ip C'/x V-'reac • 



(4.38) 



Inserting Eq. ( [4.25| ), or better, Eq. ( f4.27| ) into Eq. ( |4.38| ) yields, after integration by parts, 
a total 4-momentum loss 

4 



plost if 

M - 3 



with integrand 



z ■ z 



z ■ z 



+ 2z 



z ■ z 



The square of vr^ reads 



Viz' (i • zf + {zy {:zy - vi' {zy {z ■ z) 



(4.39) 



(4.40) 



(4.41) 



This is negative definite in the temporal gauge t = t, showing that |P|^''*| < -E", 
physically expected. 



"■lost 



as 



E. Gravitational and axionic radiation reaction 



We are going to see that the generalization of the dilaton results to the case of the 
gravitational and axionic fields is non-trivial, and leads to physically nonsensical results. 
Let us first generalize Eq. ( [4.381 ). The relations 

f/"^ dahp, = -daiy^r' K, ^6^°) + K, 5,(^7 7'^' ^b^^") , (4.42) 
dxB,^ = dkie''' daz" B,,) - ^^(e'^' 9,^^) B,, , (4.43) 
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show that, as far as their integrated effects are concerned, the gravitational and axionic 
reaction forces ( |2.1tj| ), ( |2.17| ) are equivalent, respectively, to: 



^hequiv. ^_^^./3^^^reac^ (4-44) 
jrB equiv. = 1 ^ ^reac _ ^4 45) 

It is important to note that, as in the dilaton case Eq. ( |4.38| ), these equivalent reaction forces 
are simple bilinear forms in the vertex operators and the derivatives of the fields. They can 
both be written as 



•^r'^- = ^^-^ ^-c (4.46) 



where, as in Eq. (|]34D, S denotes the source of the field A = hap or -8^/3, and where the dot 



denotes a certain symmetric bilinear form acting on symmetric or antisymmetric tensors. 
With the normalization of Eq. (|4.46| ) these bilinear forms are, respectively. 



Uo.p = Upo.: U-U = U^p U'^^ -\u^ = ' (4-47) 

K/3 = -V>a: V-V = V^pV''r (4.48) 

One can recognize here the quadratic forms defined by the residues of the gauge-fixed prop- 
agators of the h and B fields. Note that if we wish to rewrite the scalar reaction force ([4.38| ) 



in the same format ( [4. 461) as the tensor ones we have to define the dot product for scalar 
sources as 

S^-S^ = 8Ej. (4.49) 

Using this notation, and the results above on the reaction fields, it is possible to compute in 
a rather streamlined way the total 4-momentum lost under the action of the local reaction 
force: 

= / JdadrE-ld^ A,J°'^' . (4.50) 

The calculation is simple if one uses the form ( [4.27| ). Let us note that the worldsheet-scalar 
sources {S = S/^/7) for the three fields we consider (yj, h and B) satisfy 

S-S = const. , (4.51) 
S-S = 0. (4.52) 
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Indeed, if we introduce the scalarized vertex operators (with conformal dimension zero) 
U = U / ^/7, Uai3 = Ua/3/ ^/y and Va/3 = Va/3/ 1/7, it is easily seen that 

ip: U -U = 8U^ = +32, 

h: U^f3- f/°^ = Uaf3 ?7°^ - i f/2 = , (4.53) 
B : V^p- = -2 . 

The relations ( ^.51| ), ( |4.52| ) simplify very much the evaluation of P^"*^*. In particular, the 



constancy of S-S allows one to integrate by parts on i;'^, etc . . . without having to differentiate 
the S ■ S factors. By some simple manipulations, using also the consequence 

S ■ S + S ■ S = 0, (4.54) 

of Eq. (WM ), we get 



1 -2 1 12 

3 12^. 



+ ^z^{S-S)<j>\ . (4.55) 



K°'' = ^J JdadTlz,[{S-S) + iS-S) 
Here we introduced a special notation for the conformal factor (Liouville field), 

ds^ = e'^{-dT^ + da^) , e'^ = ^ = z'^ = -i" , (4.56) 
and we defined the unit time like vector and its first derivative 

^. = ^(e-^/^i.), ^^ = -|^ + ^0^>O. (4.57) 

Let us now prove the remarkable result that the contribution proportional to S'-S' in Eq. (|4.55| ) 
vanishes for all three fields when evaluated (as we are iteratively allowed to do) on a free 
string trajectory: 

" da dr z^{S ■ S) = 0. (4.58) 

Indeed, for the scalar case U = —2 and {7 = 0, while for the other fields a straightforward 
calculation gives 

U^pU --U =-a^^, (4.59) 

Va(3 V''^ = , (4.60) 

when taking into account the vanishing of terms proportional to the worldsheet derivatives of 
2;^ = —z'^ + z^". [These results have a nice geometrical interpretation linked to the Gauss- 
Codazzi relations.] Integrating by parts, we see that the contribution ( |4.58|) is proportional 
to J J da dr (□,, z^) which vanishes, again because of the free string equations of motion. 
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Finally, remembering the constancy oi S ■ S, we get the very simple result 

3 ~~8G 

where the integrand 



C = ^A^ / IdadTTT,, (4.61) 



7if, = z^ (u^ + ^^^^+z^^ (4.62) 
is easily seen to coincide with the one which appeared above, Eq. (|4.40| ), in our direct 



calculation of the dilaton reaction. Let us recall that the present calculation applies uniformly 
to all three fields if we define the dot product between dilatonic vertex operators with an 
extra factor 8, see Eq. ( [4.49| ). 



The conclusion is that the local approximation to back reaction for the three fields 
h and B leads to energy-momentum losses which are proportional to the same quantity 
J J (P (jTTfj^ with coefficients respectively given by (using Eqs. ( [4.53[ ) above) 

J = J = +^ A G a V^ (4.63) 

The result ( [4.63|) coincides with Eq. ( [4.39| ) above (for which we have verified that the overall 
sign is correct). We therefore conclude that the "local reaction approximation" ( [4.18|) yields: 
(i) a vanishing, net energy-momentum loss for the gravitational field, and (ii) the wrong sign 



(antidamping) for the axionic field. The latter result disagrees with Refs. |T^, []T5[ (see the 
Appendix) which claimed to obtain positive damping. It is for clarifying this important sign 
question that we have presented above a streamlined calculation showing that the overall 
sign can simply be read from the contraction of the vertex operators of the fields. Indeed, 
finally the physical energy-loss sign is simply determined by the easily checked (and signature 
independent) signs in Eqs. ( |4.53| ). 

F. Gauge invariance and mass-shell-only positivity 

Why is the "local back reaction approximation" giving physically unacceptable answers 
in the cases of gravitational and axionic fields but a physically acceptable one in the case of 
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the dilatonic field? The basic reason for this difference between /i^j, and on one side, 
and if on the other is the gauge invariance of the former. Indeed, a gauge symmetry (here 
h^u hf^^ + d^^i^ + du C,^, B^y — > B^y + Ay — dy A^) means that some of the components of 
h^y and B^y are not real physical excitations. This is associated with the fact that some of 
the components of h^y and B^y (namely h^i and i?oi) have kinetic terms with the wrong sign, 
i.e. that they (formally) carry negative energy. Therefore, approximating radiation damping 
is very delicate for gauge fields. A slight violation of gauge invariance by the approximation 
procedure can lead to antidamping (the literature of gravitational radiation damping is full 
of such errors, see e.g. 0]). A more precise way of seeing why the local back reaction 
approximation is dangerous in this respect is the following. 

We have proven above that an exact expression for the 4-momentum of the source lost 
to radiation is given (for ip, h^y and B^y, and more generally for any linearly coupled field) 
by an expression of the form 

= -kj d'^X J{X) ■ d, Areac(x) , (4.66) 

where J{x) is the source of A{x) 

aA{x) = -J{x), (4.67) 

and where is a positive coefficient which depends on the normalization of the kinetic terms 
of A{x) [Ank = 1/8G when using the above normalizations, the extra factor An compensating 
for our present way of writing the field equation ( [4.67|) .] The spacetime source J{x) is linked 



to our previous string distributed sources by J{x) = An J (Pa S (5^(x — z). The dot product 
in Eq. ( |4.66| ) is the symmetric bilinear form defined in Eqs. ( [4.47] ), ( [4.48| ), ( [4.49 ) above for 



the three cases h, B and ip. Introducing Fourier transforms, with the conventions, 

J{p) = [ d^xe-'P'^Jix), (4.68) 



acac(x) = ^ [aet(x) - Gadv(x)] = J ^ G,,Up) e+*^" , (4.69) 



the energy loss ( |4.66| ) reads 



PT = -kj^. GrUp) J{-P) ■ Hp) ■ (4.70) 



i2n) 
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To see the positivity properties of Pj^^^ we need to insert the exphcit expression of the Fourier 
transform of Greac- 

The Fourier decomposition of the retarded and advanced Green functions (DC = —5^) 
read 

where t] is any positive infinitesimaL Using the formula 

= P^—^in6(x -a) , (4.72) 

X — a±ir] X — a 

where P denotes the principal part, one finds 

GreM = \ \GrM - Gadv(p)] = ^^T sigU (/) ^(p^) . (4.73) 



Inserting ( [4. 731 ) into ( f4.7(J| ) one gets 

= + fcvr I ^ sign (/) Sip') J{-p) ■ J{p) 

= +kf dpp,r{p)-J{p), (4.74) 

JV+ 

where V+ denotes the positive mass shell p^ = +v^p^ and dp = {2tt)~^ d'^p/2p^ the natural 
integration measure on Vl^.. Here, we have used the reality of the source: J*{x) = J{x) =^ 
J*{p) = J{-P)- 

As in the case of Eq. ( |4.33| ) and its kin, the meaning of Eq. ( [4.74| ) is formal when 



evaluated on a (periodic) free string trajectory. However, it is, as usual, easy to convert Eq. 
( [4.74| ) in a result for the average rate of 4-momentum loss by using Fermi's golden rule: 



[(5(p° - ntu)f = ^ 5(/ - ncu) I dr . (4.75) 



2n 

One then recovers known results for the average energy radiation from periodic string motions 

The integrand in the last result has the good sign (i.e. defines a vector within the future 
directed light cone) if the dot product J*{p) ■ J{p) > 0. This is clearly the case for a scalar 
source, but for the gauge fields h^^ and B^^ one has integrands 

1 
2 



Jh,Ap) Jrip) - 7^ \ jLip)\\ (4.76) 
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and 



J*B,uip)JS'ip) (4.77) 

which are not exphcitly positive because of the wrong sign of the mixed components Joi- 
As is well known this potential problem is cured by one consequence of gauge invariance, 
namely some conservation conditions which must be satisfied by the source. In our case the 
gravitational source Jj^^^x) oc T^u{x) must satisfy d'^ J^^, = 0, while the axionic source must 
satisfy d'^ J^A^) = 0. In the Fourier domain this gives Jjiuip) = or p'^ J^uiv) = 0- These 



transversality constraints are just enough to ensure that the integrands ( [4.76| ), ( [4 .771 ) are 



positive when evaluated on the mass shell V^. What happens in the "local back reaction 
approximation" is that one replaces the Green function G'rcac(a^) by a distributional kernel 
^100(3;) with support (in x space) localized at x = 0. Its Fourier transform G'ioc(p) is no 
longer localized on the light cone = 0, and therefore the delicate compensations ensuring 
the positivity of the integrands ( 4.76|) , ( [4.77 ) do not work anymore. This explains why the 



local back reaction approximation is prone to giving unreliable expressions for the damping 
due to gauge fields. On the other hand, in the case of a scalar field the crucial source 
integrand J*{p) J{p) in Eq. ( |4.74| ) remains positive-definite even off the correct mass shell. 



This explains why, in the case of the dilatonic field, the local back reaction approximation 
might (as it was found above to do) define a physically acceptable approximation to the 
exact, non-local damping effects. 

V. IMPROVED DILATONIC REACTION AS SUBSTITUTE TO 
GRAVITATIONAL REACTION 

As the main motivation of the present study is to find a physically reasonable, and 
numerically acceptable, approximation to gravitational radiation damping, the results of 
the previous Section would seem to suggest that the local back reaction approach fails to 
provide such an approximation. However, we wish to propose a more positive interpretation. 
Indeed, both the direct verification of Section |iV D| , and the argument (in Fourier space) of 



Section |IVF| shows that the local back reaction approximation can make sense when applied 
to scalar fields. On the other hand, Damour and Vilenkin in a recent study of dilaton 
emission by cosmic strings have found that, in spite of their genuine physical differences, 
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gravitational radiation and dilatonic radiation from strings are globally rather similar. For 
the samples of cuspy or kinky loops explored in Ref. [^], the global energy losses into these 
fields turned out to be roughly proportional to each other. Even when considering in more 
detail the physically important problem of the amount of radiation from cusps, it was found 
that (despite an expected difference linked to the spin 2 transversality projection) both 
radiations were again roughly similar. 

Let us also recall that this similarity, or better brotherhood, between gravitational and 
dilatonic couplings is technically apparent in the similarity of their vertex operators (which 
are both subsumed in the form (^i, d"" z'^ da with a generic symmetric polarization tensor 
C,^nu) and is a very important element of superstring theory. This leads us to propose to use, 
after a suitable normalization, the physically acceptable local dilatonic back reaction force 
as a substitute for the gravitational radiation one. In other words, we propose to use as 
"approximation" to gravitational radiation damping a local reaction force of the form (in 
conformal gauge) 

4 



. f z ■z\ I ( z' -z 



(5.1) 



We note also that, though there are more differences between axionic and gravitational 
radiations than between the dilatonic and gravitational ones, they are still roughly similar in 
many ways (as witnessed again by the brotherhood of their vertex operators C,^^ d°- z^^ da z^ 
with now a generic asymmetric polarization tensor) so that one can hope to be able also to 
represent in an acceptable manner axionic radiation damping by a force of the type (|5.1|) 



with the replacement G/i^ GX^ and another, suitable choice of A. [Actually, due to their 
sign error, this last proposal agrees with the practical proposition made in Refs. Jl^, [p!5[] .] 
It remains to clarify the choice of A in Eq. (|5.1|) . Up to now we have implicitly assumed 
that A was constant. There are, however, several reasons for suggesting a non-constant 
A. The first reason concerns energy-momentum losses associated with cusps. To see things 
better, let us use a temporal gauge t = t and concentrate on the energy loss implied by Eq. 
( ^.1|) . One finds simply 

= -^Gfi^J JdadrA , (5.2) 

where v(o", r) = z(o", r). At a cusp v^(o", r) = 1. As v^(o", r) < 1 everywhere, near a cusp 
one will have v^(cr, r) = 1 — (aa^ + bar + cr^) + 0{{a + r)^) where the parenthesis is a 
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positive definite quadratic form. This shows that, if A is constant, the integral E^°^^ ~ 
/ J da dr {1 — v^)~^ is logarithmically divergent (as we explicitly verified on specific string 
solutions). As the real energy loss to gravitational or dilatonic radiation from (momentary) 
cusps is finite, this shows that Eq. { f).l\) overestimates the importance of back reaction due to 
cusps. In other words, if one tries to complete the equations of motion of a string by adding 



the force (|5.1| ) with A = const, this reaction force will prevent the appearance of real cusps. 



As the calculations of Ref . |Tl| , using the "exact" non-local gravitational radiation, find that 
cusps are weakened but survive, it is clear that one must somehow soften the "local" force 
( ^.1| ) if we wish to represent adequately the physics of cusps. At this point it is important 
to note that the proposal ( ^.1|) lacks worldsheet covariance, which means, on the one hand, 
that A has introduced a local coordinate length or time scale on the worldsheet, rather than 
an invariant interval and, on the other hand, that one must specify a particular time-slicing 
of the worldsheet. As the ratio between coordinate lengths and times and proper intervals is 
locally given by the square root of the conformal factor e"^ = z'"^ = (= 1 — in temporal 
gauge), it is natural to think that a better measure of the coordinate interval A to use in 
Eq. ( |5.1| ) might vary along the worldsheet because it incorporates some power of e*^. This 
might (if this power is positive) prevent the logarithmic divergence of the integral ( ^.2|) . At 



this stage, a purely phenomenological proposal is to take A in Eq. (|5.1| ) of the form 

A{a,r)=f{-zy2X, (5.3) 

where / is a dimensionless factor, rj is a. positive power, and A the wavelength of the radia- 
tively dominant mode emitted by the string. We introduced a factor two for convenience 
because, in the case of loops for which the fundamental mode is dominant, the wavelength 
is L/2 where L is the invariant length of the loop. On the other hand, if we consider a loop 
carrying mainly high-frequency excitations, or an infinite string, it is clear that A should not 
be related to the total length L, but to a length linked to the scale of the principal modes 
propagating on the string. 

Let us briefiy comment on the lack of worldsheet covariance of Eq. ( ^.1|) and on its 
consequences. Eq. emerged as a local approximation to an integral which had the 

same formal expression in all conformal gauges. The "local back reaction approximation" 
procedure has, among other things, violated the formal symmetry between r and a on the 
worldsheet. From the formal point of view this loss of symmetry is certainly unpleasant and 
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it would be nicer to be able to write a local force density which respects the symmetry of the 
worldsheet conformal gauges, and does as well as Eq. ( |5.lD in entailing a positive energy loss 
quantitatively comparable to the result ( |5.2| ) (which will be seen below to be an adequate 
representation of the actual energy loss). We failed to find such a covariant local force density. 
This is why we propose to use (|5.1|) , despite its formal imperfections, as a substitute to the 
exact, non-local gravitational radiation damping. From this point of view, the asymmetry 
between r and a in Eq. can be interpreted as a sign that the purely local expression 
tries its best to incorporate the, in reality, global damping effects by selecting special 
time-slicings of the worldsheet (r = const, lines, and their orthogonal trajectories). A natural 



physical choice of special time-slices (necessary to define properly the meaning of (|5.1| )) is 
to consider the spatial sections associated to the (instantaneous) center of mass-frame of the 
string. [Note that the numerical calculations below of the energy loss ( |5.2| ) are performed in 
the string center of mass frame.] For a free Nambu-Goto string in flat space, this definition 
is compatible with using a worldsheet gauge which is both conformal and temporal (i.e. 
r oc P'^^^'^^z^). Therefore, in such a case, the (orthogonal) worldsheet vector field d/dr is 
well defined (both in direction and in normalization), which means that JF^, Eq. ( |5.1| ) is 
well defined, on the worldsheet, as a spacetime vector locally orthogonal to the worldsheet. 
We shall admit that the definition of JF^ can be smoothly extended to the case where the 
(Nambu) string moves in a curved background spacetime (say a Friedmann universe). When 
working in the approximation of a flat background the expression ( |5.1| ) can be used directly 
as right-hand side of the standard, flat-space, conformal gauge string equations of motion: 
'^^lv {'^^ — z"^) = T^. [Note that we use here a flat-space worldsheet gauge, Eq. (|2.7|) with 
Qiiv Vf^u-] Note finally that the actual numerical simulations of a string network introduce 
a particular time slicing and one might also decide (for pragmatic reasons) to use it to define 
the T = const, slices of Eq. ( ^.1| ) (i.e. to neglect the Lorentz-transformation effects associated 
to the center of mass-motion of the strings). 

A first check of the physical consistency of the proposal ( ^.3| ) consists in verifying that, 
despite the nonconstancy of A, the integrated energy loss ( ^.21) will be positive for all possible 
loop trajectories. Integrating by parts Eq. (p.2|) one finds 



= lG/iV(2A) / I dadT{l-vy 



V ■ V , . (v ■ v)^ 



+ 2(1-7^) 



(5.4) 



1 - V2 ' ' " (1 - V2)2_ 

This is manifestly positive (and finite) as long as < < 1. Assuming this to be the case. 
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the question is then: Are there values of / and r] (after having decided on a precise definition 
of A) such that the corresponding damping force ( |5.1| ) gives a reasonably accurate description 
of the "exact" effects of energy loss to gravitational radiation? We did not try to answer this 
question in full generality. For simplicity, we fixed the power rj to the value f] = \ (which 
seems intuitively preferred as it evokes a Lorentz contraction factor arising because we look 
at string elements "moving" with relativistic speeds). Then we compared the energy loss 
due to ( p.lD to the energy radiated in gravitational waves as (computed (using Eq. ( [4.74[ )) 
in the literature (both energy losses being evaluated in the rest frame of a free string). As a 
sample of loop trajectories we consider Burden loops 

1 



z(r,(T) = - [a{u)+h{v)] 



(5.5) 



L 

L 
2^ 



— cos(m m) 6*3 H sin(m u) e\ 

m m 

— cos(n f ) 63 sin(n v) e[ 

n n 



where 



27r 



u 



2ti . 



COS 6*1 + sin 1p 6*2 . 



(5.6) 
(5.7) 

(5.8) 



This family of solutions depends on the overall scale L, which is the total invariant length 
of the loop (M = fiL), on two integers m and n, and on the angle ip. Our parameter ip 
coincides with the angle ip in |^^, denoted in f^. The actual oscillation period of the loop 
is T = L/{2mn) which leads us to choosing 2A = 2T = L/mn in Eq. ( ^.31 ). With this choice 
we computed the energy loss (|5.2| ). The calculation is simplified by noting, on the one hand. 



that, for this family of loops, v ■ v = — |^ m^+n^ ^ ^^^q other hand that the 

worldsheet integral in ( ^.2|) can be rewritten in terms of an average over linear combinations 
of the two angles 27rm(r — a)/ L and 27™(r + a)/L. This yields simply for the average rate 
of energy loss 



r = 



E 



lost 



/ 



2mn 



■7 



where 
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FIG. 1. Dimensionless energy loss rate Burden loops with (m, n) = (1, 1) and (1, 3). 



with 



1 — - (1 + cosip) cosx — - (1 — cosip) cosy 



(5.11) 



We plot in Fig. |I| Tm,n as a function of the angle ip, for the nominal value / = 1 and for the 
two cases (m, n) = (1, 1), (m, n) = (1, 3). [As said above there is a simple scaling law for the 
dependence on m and n.] If one compares this Figure with the figures published in |33], [0 
(Fig. 7.6, p. 205 there) one sees that they give a roughly adequate numerical representation 
of energy losses to gravitational radiation if 



/^0.8. 



(5.12) 



The fact that our present "best fit" value of the factor / leads to values of A which are 
numerically comparable to L (when {m,n) = (1, 1)) rather than to a smaller fraction of L 
should not be considered as physically incompatible with the idea of using a local approx- 
imation to back reaction. Indeed, on the other hand, the rough justification of the local 
approximation given in Section |TV| suggested A ~ n/u ~ A/2, i.e. something like L/4, and. 
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on the other hand, numerical computations show that the energy lost to dilaton waves (with 
coupling a = 1) is smaller than that lost to gravitational waves by a factor of order 3 or 
so (part of which is simply due to the fact that there are two independent tensor modes 
against one scalar mode). Therefore, as we use A only as an effective parameter to model 
gravitational damping it is normal to end up with an increased value of A/L. 

Clearly, more work would be needed to confirm that the modified local dilaton reaction 
( p.lD can be used as a phenomenological representation of gravitational reaction. Our main 
purpose here was to clarify the crucial sign problems associated to gauge fields, and to 
give a first bit of evidence indicating that Eq. ( ^.1|) deserves seriously to be considered 
as an interesting candidate for mimicking, in a computationally non-intensive way, the back 
reaction of gravitational radiation. We are aware that several important issues will need to be 
further studied before being able to use Eq. ( |5.1| ) in a network simulation. Some numerically 
adequate definition of A will have to be provided beyond a case by case definition, which 



in the case of long loops decorated by a regular array of kinks, as in Ref. |Tl|, would be 
something like 2 A ~ L/N where N is the total number of kinks. We note in this respect 
that a Burden loop with m = 1 and n ^ 1 provides a simple model of a long, circular 
loop decorated by a travelling pattern of small transverse oscillations. However, the local 
approximation (|5.1| ) cannot expected to be accurate in this case, because the radiation from 
purely left-moving or right-moving modes is known to be suppressed 0. This suppression 
is not expected to hold in the more physical generic case where the transverse oscillations 
move both ways. The accuracy of the local approximation (|5.1| ) should therefore be tested 
only in such more generic cases. 

The explicit expression ( |5.1|) must be rewritten in the temporal, but not necessarily con- 
formal, worldsheet gauges used in numerical simulations, and the higher time derivatives in 
jF^j must be eliminated by using (as is standard in electrodynamics and gravitodynamics 



1 3^) the lowest-order equations of motion. [These last two issues have already been treated 
in Refs. [^, |]15[-] Finally, we did not try to explore whether 77 = 1/2 is the phenomenolog- 
ically preferred value. To study this point one should carefully compare the effects of (|0| ) 
on the weakening of cusps and kinks with the results based on the exact, non-local reaction 
force [|Tl|]. [The facts that the curves in Fig. ^ are flatter than the corresponding figures in 
[0] suggest that a smaller value of r] might give a better fit.] 
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VI. CONCLUSIONS 



In this paper we studied the problem of the radiation reaction on cosmic strings caused 
by the emission of gravitational, dilatonic and axionic fields. We assume the absence of 
external fields. We use a straightforward perturbative approach and work only to first order 
in G. Our main results are the following. 



Using the results of Refs. [|T^, [|1^] for the renormalization of the string tension /i, 
we write down the explicit form, at linear order in G, of the renormalized equations 
of motion of a string interacting with its own (linearized) gravitational, dilatonic and 
axionic fields. [Within our framework, we verified the on shell finiteness of the bare 
equations of motion, which is equivalent to their renormalizability.] 

We have extended a well-known result of Dirac by proving for general linearized 
fields that, in the decomposition ( [4.14| ) of the renormalized self-force, only the time- 
antisymmetric contribution J^^^^^ = jF^(Areac), where A^-cacix) is the half-retarded mi- 
nus half-advanced field, contributes, after integration over time, to the overall damping 
of the source. [This result had been assumed without proof in previous work on the 
topic] The "reactive" self-force J^J^^^ is manifestly finite (and independent of the 
renormalization length scale Ar), and is non-local. 

We have critically examined the proposal of Battye and Shellard [l^, [|l5l (based 



on an analogy with the Abraham-Lorentz-Dirac treatment of self-interacting point 
charges) to approximate the non-local integral ( [4.15|) entering the reactive self-force 



jrveac ^y. local expressiou ( [4.18|) . For this purpose we found very convenient to 



use dimensional continuation, a well known technique in quantum field theory. We 
found that the local back reaction approximation gives antidamping for the axionic 
field, and a vanishing net energy-momentum loss for the gravitational one. We argued 
that the ultimate origin of these physically unacceptable results come from trying to 
apply the local back reaction approximation to gauge fields. The non-positivity of 
the local approximation to the damping comes from combining the modification of the 
field Green functions implicit in the local back reaction method, with the delicate sign 
compensations ensured, on shell only, by the transversality constraints of the sources 
of gauge fields. 
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By contrast, we find that the local approximation to the dilatonic reaction force has 
the correct sign for describing a radiation damping. In the case of a non-gauge field 
such as the scalar dilaton there are no delicate sign compensations taking place, and 
the coarse approximation of the field Green function, implicit in the local back reaction 
method, can (and does) lead to physically acceptable results. 

Taking into account the known similarity between the gravitational and dilatonic ra- 
diations (e.g. P), we propose to use as effective substitute to the exact (non-local) 



gravitational radiation damping the "dilaton-like" local reaction force ( |5.1| ), with a 
suitably "redshifted" effective length A, Eq. (^.3|). This force is to be used in the 
right-hand side of the standard, fiat-space conformal-world-sheet-gauge string equa- 
tions of motion, with r- slicing linked, say, to the global center-of-mass frame of the 
string. The numerical calculations exhibited in Fig. |I] give some evidence indicat- 
ing that Eq. ( ^.1|) deserves seriously to be considered as an interesting candidate for 
phenomenologically approximating, in a computationally non-intensive way, the back 
reaction of gravitational radiation. [We recall that the exact, non-local approach to 
gravitational back reaction, defined by Eq. ( [4.14[ ), is numerically so demanding that 



there is little prospect to implementing it in full string-network calculations.] More 
work is needed (e.g. by comparing the dynamical evolution of a representative sample 



of cosmic string loops under the exact renormalized self-force ( [4.14| ) and our proposed 
( p. ID ) to confirm that our proposed substitute ( p.lD is a phenomenologically acceptable 
representation of gravitational reaction (or of the combined dilatonic-gravitational re- 
action, as string theory suggests that the dilaton is a model-independent partner of 
the Einstein graviton). 

It will be interesting to see what are the consequences of considering the effective reaction 
force, Eq. ( |5.1| ), in full-scale network simulations (done for several different values of Gfi) 
of gravitational radiation. Until such simulations (keeping track of the damping of small 
scale structure on long strings) are performed, one will not be able to give any precise pre- 
diction for the amount and spectrum of stochastic gravitational waves that the forthcoming 
LIGO/ VIRGO network of interferometric detectors, possibly completed by cryogenic bar 
detectors, might observe. 
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APPENDIX A: 



In this Appendix we will give some details on the derivation of Eqs. ( [4.24| ), ( [4. 251 ) using 
dimensional continuation. 

A nice feature of analytic continuation is that it allows one to work "as if" many sin- 
gular terms were regular. For instance, the factors and that appear 
in Eqs. ( [4.22|) , ( [4.23| ) blow up on the light cone (fi^ = 0) when n = 4. However, if we 
take the real part of e = 4 — n large enough (even so large as corresponding to negative 
values for Re(?T,)), these fi-dependent factors become finite, and actually vanishing, on the 
light cone. This remark allows one to deal efficiently with the f2-dependent factors ap- 
pearing in Eqs. ( [4.22|) , ([4.23|) . We are here interested in the contributions to Aj.(,t{z) and 
dfj, Aj-etiz) coming from a small neighbourhood z' = z{t',(j') oi z = z{T,a) on the world- 
sheet. Let us, for simplicity, denote uj = VP'. We first remark that when (t', a') {ti<^)i 
ijj = —{z{r, a) — z{t', a')Y admits an expansion in powers of r' — r and a' — cr of the form 



\a — a 



UJ = UJ2 + ^3 + + ■ ■ ■ , 

^1, and 



(Al) 



with UJ2 = —z [{t' — t) 

uj^ = ((r' - rf + (r' - rf {a' - a) + (r' - r) (a' - af + (a' - af) , etc . . . (A2) 



Then we can formally expand the f2-dependent factors of Eqs. ( [4.22| ), ( ^.23| ) in powers of 
t' — T and cr' — cr as follows 



T [tu" d{uo)] 



e{uJ2) + ^(cUs) (CU3 + CU4 + . . .) + 6\UJ2) {uJ^i + . . .f + . . . 



X 



Here and below, the symbol T will be used to denote a (formal) Taylor expansion of any 
quantity following it. This expansion is valid (at any finite order) when Re(a) is large 
enough, and is therefore valid (by analytic continuation) in our case where a = {2 — n)/2 or 
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—n/2. A technically very useful aspect of the above expansion is that all the terms containing 
5{ijJ2) or its derivatives give vanishing contributions (because uj2~^ 5^^\ijJ2) vanishes if Re(a) 
is large enough, so that, by analytic continuation, uo2~^ 5'^^\uj2) = for all values of a). The 
net effect is that the contribution coming from a small string segment —^<{a' — a)<^ 
around a (with A being much smaller that the local radius of curvature of the worldsheet) 
can be simply (and correctly) written as the following expansion: 

ro-+A/2 



[A,,,{z)f = e da' dT'T.uj—9{uj) 

= e / da' dr'T (T.uj—] e{uj2) 

J(7-A/2 J~oo ^ ' 



\/2 

:/ da' dr'T (j:cu—) 

Ja-A/2 J-oo ^ ^ 



/cr-A/2 

Here, we have introduced an arbitrary upper limit ri, submitted only to the constraint 
Tret < Ti < Tadv (for iustaucc Ti could bc t), and which replaces the missing theta function 
6{z^ — z'^) by selecting the retarded portion of the other theta function 9{uj). As above, the 
symbol T denotes a formal Taylor expansion. The expansion T(Sco'") is simply obtained by 
multiplying the expansion (|A1|) of uj with that of S(t', a'), namely 

T [S(r', a')] = S(r, a) + (r' - r) E + (a' - a) S' + . . . (A3) 

Similarly we have 

A,,,{z)f = 2e / da' / dr' T (S uj-^'^) , (A4) 

J0--A/2 J-oo 



rcr+A/2 _ fT-\A(j 
la-A/2 

as well as corresponding expressions for the advanced fields 

[^adv(^)]^ = e / da' dr'T (Sc.— , (A5) 

/•cr+A/2 /■+00 

^adv(^)]^ = 2e / da' / dr' T (s u"'/') . (A6) 

Jcr-A/2 JT+\Aa\ ^ ' 

As a check, we first computed the ultraviolet divergent contributions to A^^'^{z) and A^^t^,^)- 
We find 

^ret(^) = -4 f-) 2S, (A7) 



(A8) 
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As it should, Eq. ( [A8D yields exactly the same divergences as we found in Sec. |T| by 
introducing a cut-ofF 6 in the a' integration in four dimensions. More precisely, Eq. ( |AS| ) 
coincides with Eq. ( |3.11| ) if we change 2/e — log 1/6. Let us note that, in the present 
approach, the renormalization scale Ar would enter by being introduced as a dimension- 
preserving factor in the dimensionful coupling constants, like Newton's constant G, say 

Our main interest is to compute the "local approximations" to the reaction field 

1 



(v4rct(x) - v4adv(a;)) 



(A9) 



and its derivatives. Dimensional continuation gives an efficient tool for computing these. 
Indeed, combining the previous expansions we can write 

A..r(x) = -e I da' I dT'e{^i)Tu'^.^..A^{^^)^) , (AlO) 



A/2 I-+00 

da' / dT'e{nl)T(^,,^,)^,, 

Jt+\A<t\ 

A/2 



da' \ dT'Q(Sll)T(^,,_,^^,, hl^^iSl^yA , (All) 

-A/2 Jt+|Ao-| ^ ' 

where T(^'_T-)odd denotes the part of the Taylor expansion which is odd in r' — r. Moreover, 
as we know in advance (and easily check) that the a'-integrands in Eqs. ([A10|) and ( |A11 



are regular at a' = 0, we can very simply write the result of the local approximation (|4.18| ) 
(with a corresponding definition for A^°^'^^{z)) by replacing a' = a in the integrands of Eqs. 



[^reac(^)] 
[f^/j ^reac(^)]'°'^'^^ - 



-eA 
-2eA 



-^(r'-r) odd 



-'(r'-r) odd 



(A12) 
(A13) 



Here, T(jK^r)oAA denotes the operation of replacing a' by a and keeping only the odd terms 
in the remaining Taylor expansion in r' — r. This simplifies very much the computation of 
the reactive terms (making it only a slight generalization of the well known point-particle 
results, as given for a general source in, e.g. ||3ll]). Indeed, inserting the following expansions 

\{r'-rfz,, (AM) 

1 



l]^(r',a)~ 



r - r z„ 



r?z,, 



S(r',a) ~ S(r,(T) + (r' - r) S 



-' -tYT.^ 



~2 / / 



Qr^T' , a) ~ —z" Kj' 
1 



1-f 






- (r - r) 











z ■ z 



z" 



r' - r)2 



z ■ z 



(A15) 



(A16) 
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in Eqs. (|A12|) , (|A13| ) we get our main results 



local 



[^reac(^)] 



A 



A 



S - S 
1 



z ■ z 



y2\2 



+ S 5^ + S i:^ - 4S z^ 



z ■ z 



2S^„ 



z ■ z 



-S z„ 



z ■ z 



3 ^ 



z ■ z 



6S z,, 



z ■ z 



(A17) 



(A18) 



These results were also obtained (as a check) from Eqs. (|A10|) , ( All ) without using in 
advance the simplification of putting a' = a in the integrand. 

We have also performed a direct check on these final expressions by comparing them to 
the well known point-particle case 



30 



3l|. Indeed, we have seen above that Aj-eac' 



and d^J^°^^l could be thought of as being generated by the effective source YT^-{t' ,a'^ = 
6{cr' — cr) A S(t', a), i.e. a source along the world-line C^j, defined by a' = a. For any given 
value of a, by transforming the coordinate time r' into the proper time s = J e'^^'^ dr' along Cfj 
and by renormalizing in a suitable way the source A S(r', cr) = e^^^^ gl^g'^ ( go that the stringy 
spacetime source / d'^a'Il'^^'{T' ,a') 6'^{x — z{a')) transforms into the standard point-particle 
source / ds S{s) 6'^{x — z{s))), we recovered from Eqs. ( [A17| ), ( |A18| ) known point-particle 
results ||31|. This check is powerful enough to verify the correctness of all the coefficients in 



Eqs. (IATtD , (|AT8 



m 



In order to compare directly our expressions with what derived by Battye and Shellard 
n|, |15|, let us write Eq. ( [A18| ) for the axion field. We get 



AGXA 

(^2\2 



3 



z ■ z 

"1^ 



z ■ z 



(A19) 



where i^t^H = /^A^ti^ _|_ ^fivx _|_ T^i^A^t^ Note that, when identifying the basic contravariant 
tensors z^ and V^^ , the tensor H^^^'^ (and the force density J^^) must be identical in our 
conventions and in the ones of Refs. |]14| . 



(who use the opposite signature). However, our 
result Eq. ( A19 ) differs, after the substitution G\ /a/8, in many terms from the second 
Eq. (31) of Ref. |jl5|. Whatever be the corrections we could think of doing on the second 
term in their Eq. (31) (which is dimensionally wrong, probably by a copying error leading 
to a forgotten overdot on one of the two terms), we saw no way of reconciling their result 
with ours (even after expanding explicitly V^'^ = z^ z^' — z^ z^'\ 
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